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FOREWORD 


The research in this document was intended to contribute to determining 
the feasibility of predicting the conditions imder which wind/ turbulence environments 
hazardous to aviation operations exist. Aeronautical Research Associates of Princeton 
has developed a computer model for solving for the velocity, temperature, and turbu- 
lence distributions in the atmospheric boundary layer. The model is based on using 
invariant modeling for clousre of the dynamic equations of the ensemble-averaged, 
single-point, second-order correlations of the fluctuating velocities and temperature. 

The model equations are reviewed, and the dimensionless parameters investigated 
include surface roughness, geostrophic wind, latitude, surface heating, baroclinicity, 
radiation flux divergence, and geostrophic accelerations. Results of some variations 
in each of these variables and of some indication of probable maximum wind shear 
conditions at airports are presented. From the National Transportation Safety Board 
aircraft accident data bank, a set of accidents for which wind shear and/or turbulence 
appeared to be a contributing factor was investigated to select cases which can be 
used to evaluate the feasibility of predicting when conditions will be hazardous to air- 
craft operations. It is believed that results from this two-part study can have signifi- 
cant implications relative to the aeronautical safety aspects during landing and takeoff. 

This research was conducted by the Aeronautical Research Associates of 
Princeton for the National Aeronautics and Space Administration, George C. Marshall 
Space Flight Center, Huntsville, Alabama, under the technical direction of Mr. Dennis 
W. Camp and Mrs. Margaret B. Alexander of the Space Sciences Laboratory. The 
support for this research was provided by Mr. John Enders of the Aeronautical Operating 
Division, Office of Advanced Research and Technology, NASA Headquarters. 
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WIND SHEAR AND TURBULENCE AROUND AIRPORTS 


W. S. Lewellen and Guy G. Williamson 


SUMMARY 


In June of 1975 j we initiated a research program with the 
ultimate goal of determining the feasibility of predicting the 
conditions under which wind/turbulence environments hazardous to 
aviation operations exist. This report on the first six months 
of this effort is divided into two parts. The first part des- 
cribes the computer model used for solving for the velocity, 
temperature, and turbulence distributions in the atmospheric 
boundary layer and gives the results of a parametric analysis to 
determine the expected range of wind shear and turbulence to be 
encountered in the vicinity of airports. The second part des- 
cribes the delineation of an ensemble of aircraft accidents in 
which low-level wind shear and/or turbulence appeared to have 
been causative factors. This set of accidents, encompassing a 
wide range of meteorological conditions, should prove useful in 
developing techniques for reconstructing hazardous wind environ- 
ments for aircraft safety investigation purposes. 

Also during this six months, an investigation of the wind 
shear and turbulence conditions existing at Kennedy International 
Airport on the afternoon of 24 June 1975 was made. Results of 
this study are detailed in a separate report. It showed that a 
wind shear of the order of 15 m/sec over a 100 m altitude change 
was quite consistent with the meteorological conditions at the 
time of the Eastern Airlines Plight 66 crash. 

PART 1. PARAMETRIC ANALYSIS OP TYPICAL PROPILES 
IN THE PLANETARY BOUNDARY LAYER 

INTRODUCTION 


Over the last few years. Aeronautical Research Associates of 
Princeton, Inc. (A.R.A.P.) has developed a computer model for 
solving for the velocity, temperature, and turbulence distribu- 
tions in the atmospheric boundary layer (refs. 1-5). The model 
is based on using invariant modeling for closure of the dynamic 
equations of the ensemble-averaged, single point, second-order 
correlations of the fluctuating velocities and temperature. 


Details of our atmospheric boundary layer model, as developed 
under EPA sponsorship for the purpose of predicting the dispersal 
of atmospheric pollutants, are given In reference 5. The basic 
assumption Is that the dependence of the third-order correlations 
of the velocity fluctuations on the second-order correlations, the 
mean flow properties, and their derivatives Is Invariant with 
respect to changes in flow geometry. This permits data from 
relatively simple flow experiments to be used In evaluating the 
necessary coefficients in the modeled terms. 

Verification comparisons between model predictions and a 
variety of experimental observations are presented In references 
5-8. The prediction of wind shear and vertical velocity fluctu- 
ations has been tested against observations in the surface layer, 
and agreement Is very good (ref. 7). Also, predictions of the 
vertical velocity fluctuations In the afternoon unstable mixed 
layer agree very well (ref. 8) with laboratory simulations of this 
layer. These verifications have convinced us that the model is 
fundamentally correct. However, specific verification of wind 
shear and turbulence prediction capabilities In the vicinity of 
airports should be tested by comparison with experimental observa- 
tions at a number of airports In order to permit us to set more 
precise confidence levels. This process might lead to some 
adjustments In the model, but we believe the present estimates 
represent the best currently available. 

Herein we will first review the model equations and then 
formulate the governing dimensionless parameters to be studied. 

The dimensional variables to be Investigated are surface rough- 
ness, geo.strophlc wind, latitude, surface heating, barocllnicity , 
radiation flux divergence, and geostrophlc accelerations. Results 
of some variation In each of these variables will be presented. 
Although In the confines of this brief study It has not been 
possible to look at all possible combinations of these variables*, 
it Is possible to give some Indication of probable maximum wind 
shear conditions at airports. 

As a part of this parameter study, we have Investigated 
conditions approximating those existing at Kennedy International 
Airport on the afternoon of 24 June 1975 when Eastern Airlines 
Plight 66 crashed. The results of that Investigation have been 
detailed In a separate report (ref. 9). 
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NOMENCLATURE 


f 

G 

g 

k 

L 

M 

P 

q 

Q 


R1 

Rl{ ) 


Ro 

t 

U 

Ug 

u 


V 


V 


g 


V 


W 

w 

w* 

z 


Coriolis parameter = 2R sin <p 

parameter measuring geostrophlc acceleration, eq. (19) 
gravitational acceleration 

von Karman's constant = -u*/(z 3u/9z) In the neutral 
surface layer 

Monln-Obukhov length = - 0^u*®/kg(^)^ 

parameter measuring the magnitude of barocllnlclty , 
eq. (l8) 

mean pressure 

square root of twice the turbulent kinetic energy 

temperature source strength denoting radiation flux 
divergence or temperature advectlon 

Richardson number 

bulk Richardson number based on the temperature and 
velocity difference between the altitude Indicated In 
( ) and the surface 

Rossby number based on geostrophlc wind conditions 

= Ug/fZo 

time 

mean velocity In the direction of the geostrophlc wind 
at the top of the boundary layer 

U component of the geostrophlc wind velocity 

fluctuating velocity In the direction of the geo- 
strophlc wind at the top of the boundary layer 

surface shear stress velocity 

mean velocity perpendicular to the geostrophlc wind 
direction at the top of the boundary layer 

V component of the geostrophlc wind 

fluctuating velocity perpendicular to the geostrophlc 
wind direction at the top of the boundary layer 

mean vertical velocity 

fluctuating vertical velocity 

free convective velocity scale = [gZj_(w6 )o/0o^ ^ 

vertical coordinate 

temperature Inversion height above an unstable layer 
surface roughness height 


3 


e 

0 

0 

0o 

A 

y, Uo» 

4 > 


parameter defining the direction of the geostrophlc 
wind shear, eq. (l8) 

dissipation function 

mean potential temperature perturbation 
fluctuating potential temperature 

absolute temperature about which the perturbation is 
taken 

model macroscale length 

stability parameters defined in eqs. (l6) and (17) 
latitude 

angular velocity of the coordinate system 


Superscript : 

denotes ensemble average 
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REVIEW OF MODEL 


For homogeneous boundary conditions at the surface and at the 
top of the boundary layer, the equations of motion for the mean 
horizontal velocity and perturbation potential temperature within 
the boundary layer may be written for high Reynolds number as 


|« = f (V . Vg) - ™ 

(1) 

|V = _f(u . Ug) . 

(2) 

If = - S 

(3) 


Even with the energy source term Q prescribed, the solution to 
this set of equations Is not determined until some functional 
dependences of turbulent transport terms uw , vw , and w0 are 
determined. Classical methods of modeling the planetary boundary 
layer (ref. 10) Involve the assumption of some type of eddy viscos- 
ity to determine turbulent transport. Our approach Is to use the 
Reynolds stress equations to directly compute the ensemble-averaged 
turbulent transport terms with Invariant modeling techniques applied 
to close the system of equations. These equations with their 
modeled terms may be written as follows (ref. 5): 

the six equations for the Independent Reynolds stresses 


= 4s^(uv sin <() - uw cos cj) ) - 2uw + 0 


3 

3z 


Aq 


3uu 

3z 



(4) 


9vv 

3t 


= - 4J^uv sin 4> - 2vw 


3z 


+ 0.3 


3z 


Aq 


3vv 

3z 


-K 


W - 


3 J 


12A 


(5) 


I 


3ww 

3t 


= 4nuw cos (f) + w6 + 0.3 


3z 



3ww 

3z 


a 

A 




( 6 ) 


3uv 

3t 


2n(vv sin (j> 


vw cos (f) 


uu sin <})) 



vw 


3z 


+ 0.3 


3z 



^ uv (7) 


3uw 

3t 


2^^(vw sin cp - ww cos (|) + uu cos cj)) 


3U 

3 


o 


+ 





( 8 ) 


3vw 

3t 


2J^(-uw sin (j) + uv cos (J>) 


o 


+ 0.3 


3z 




( 9 ) 


the three equations for the components of the heat flux vector 


= 2J^(v0 sin (j) - w6 cos 4>) “ 


+ 0.3 


3z 


Aq 


3u9 

3z 
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( 10 ) 


ave 

at 


-2fiu9 sin 4> - vw - w6 -5— 

dz dz 


+ 0.3 


a 

az 



0.75q 

A 


v9 


( 11 ) 


aw9 

at 


2J^u9 cos (p 





+ 



aw9 

az 


0.75q 

A 


w9 


( 12 ) 


the equation for the temperature variance 


a9 ae . n o 9 

_ = _2w9 ^ + 0.3 ^ 


Aq 


a9 


2 -1 


az . , 

j J 




(13) 


and the equation for the turbulent macroscale 


aA 

at 


0.35 


A 


— au 


+ vw 


az J 


+ 0.075q + 0.3 




az j 


0.375 /9qA\^ , 0.8A gw9 
q Uz J q2 0^ 


(14) 


These fourteen equations are solved simultaneously for the 
turbulence field. The detailed accounting of the modeled terms 
and the choice for the coefficients are given In references 3 and 
5 . Similar modeled equations are being used by Wyngaard et al. 
(refs. 11 - 13 ) and Mellor and Yamada (refs. 14-15) to compute 
atmospheric boundary layer flow. The primary differences In these 
models are that the model given In eqs. (4) through (l4) has been 
developed from comparison with laboratory data to predict a wider 
class of f lows J rather than being specialized to apply only to the 
atmospheric boundary layer flow. There are two reasons for this. 
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First, laboratory flows can. In general, be more carefully speci- 
fied and measured. Second, the wider the class of flows for which 
the model is tested, the more likely the model is to be valid for 
an untested flow situation. 

The general boundary conditions applied to these equations at 
the top of the boundary layer are that as z ^ ~ 


U = Ug 
V = 0 


geostrophic wind 


9z 


constant (difference between the adiabatic lapse 
rate and the ambient lapse rate) 


u. u. = U.6 = 60 = A = 0 

X J J. 


Rather than compute all the way to the surface with these equations, 
which would, in general, require computing the flow around various 
roughness elements, the surface boundary conditions are applied at 
Zq , the effective aerodynamic roughness height. At z = Zq , 

U = V = 0 


surface 


w6 = 


surface ^ 


and the surface layer, Monin-Obukhov functions hold for the turbu- 
lent correlations. 

Within the confines of this report, the model will be 
restricted to a one-dimensional, unsteady model with the above 
boundary conditions a function only of time. This sets the major 
limitation on the application of the model. As long as the surface 
conditions are reasonably homogeneous, the model may be applied to 
any location just by changing the boundary conditions. However, 
strong spatial inhomogeneities in surface roughness or heat flux 
will Invalidate this one-dimensional model. An extended two- 
dimensional version of this model is now being used to make calcu- 
lations along a shoreline for the Navy. 
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GOVERNING DIMENSIONLESS PARAMETERS 


The flow is now completely specified when initial conditions 
on all the variables are specified. The dimensionless parameters 
governing the flow may be obtained by normalizing the equations 
and boundary conditions. For a characteristic velocity, let us 
choose Ug . Then, If we normalize all lengths by Zq and time 
by f“S eqs . (1) and (2) Introduce one parameter, a Rossby number, 
Ro = Ug/zQf . When the flow Is neutral' (i.e., 0 = constant 

everywhere), then the only other parameter Introduced by the 
remaining equations, eqs. (4) - (9) and (l4). Is cot (|) . We will 
explore the variation resulting from changes In the Rossby number 
but will not explore the variation caused by changing cot (p 
alone. Shir (ref. l6) has reported that varying cot <|) has a 
weak influence on the correlation distributions. We believe this 
to be the case but have not made a specific study of the sensi- 
tivity to this parameter. 

• 

When the potential temperature Is allowed to vary but Qj, 
remains zero, one or more stability parameters must be Introduced. 
If the temperature Is normalized by A0, not yet specified, then 
the equations introduce only one more parameter, a bulk Richardson 
number, 

Rij, = gA0Zj,/e^u| (15) 

For stable conditions, A0 may be taken as the temperature differ- 
ence across the boundary layer, while for unstable conditions. It 
would appear reasonable to take It as the height of the Inversion 
layer times (dd/dz)^ . We will choose a bulk Richardson number 
based on the temperature and velocity differences between the 
surface and 10 meters. This appears adequate to characterize a 
quasi-steady turbulence field under stable conditions. However, 
unstable conditions appear to require two stability parameters 
even under quasi-steady conditions. To the bulk Richardson 
number, we add the ratio of the inversion height, , to Ug/f . 

For this purpose, Zj_ Is defined as the altitude near the top of 
the boundary layer where the temperature variance reaches a local 
maximum. 

This combination of stability parameters is not unique. 
Zilltinkevlch (ref. 17) uses two parameters 

g(wF) k^ 

Pq = 2 (l6a) 

* o 

and 


9 


y 


(i6b) 





Further, he argues that in relatively strongly unstable boundary 
layers, the flow becomes a function of y only. Some of our 
results will be presented in terms of these variables. 


When Q is permitted to represent some radiation flux diver- 
gence Qj. , then another governing parameter exists. One possible 
grouping for this parameter is to replace w0 in Uq with Q^Zj^ 
to give 


.l^r 




( 17 ) 


This should be a rough measure of the stability Influences on the 
turbulence resulting from radiation flux divergence. 

The remaining parameters Involve changes in the boundary 
conditions. To investigate geostrophlc wind shear, otherwise known 
as thermal wind, or barocllniclty , a variation in Ug and Vg as 
a function of z is Imposed as a boundary condltion^on the flow. 
This requires two parameters to specify the magnitude and direction 
of the geostrophlc wind shear. Following Arya and Wyngaard (ref. 
l8), these will be taken as 


1/2 z 


® @z=z^ 

_1 9 ^ 0./ 9 2 

^0 ^ 8U^/9z 

g 



(18) 


However, we have replaced u* by Ug in forming M . 

The one remaining parameter to be investigated in this report 
is acceleration in the geostrophlc wind. The most natural para- 
meter with which to compare is the Coriolis frequency f , i.e.. 


= U /|U If 
' max ' 


g 


( 19 ) 
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The approximate range of values for the 8 parameters which 
may be expected in the vicinity of airports is given in Table I. 
Sample results of the model predictions for variations in these 
parameters are given in the next section. 

Table I. 

Range of Expected Values of Planetary Boundary Layer Parameters 


Ro 

O 

H 

- 

lo8 

^^(10) 

-0.1 

- 

+ 0.1 

^0 

-500 

- 

+100 

U 

-500 

- 

+100 

Ur 

-100 

- 

+100 

M 

0 

- 

10 

3 

0 

— 

2tt 


- 10 

— 

+ 10 
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RESULTS OP PARAMETER VARIATION 


Rossby Number Variation 

Perhaps the simplest variation to unravel Is that due to 
Rossby number under neutral, barotroplc, steady-state conditions. 

The Rossby number variation may be due to changes In surface rough- 
ness, geostrophlc wind, and/or latitude. Figures 1 and 2 show the 
mean wind profiles normalized by the geostrophlc wind as a function 
zf/Ug for three different values of Ro . The wind component 
parallel to the geostrophlc wind has a small overshoot which Is 
slightly amplified as Ro Increases. The dimensionless height at 
which It occurs Is also reduced. The crosswind perpendicular to 
the geostrophlc wind also Increases and occurs at a lower altitude 
as Ro Increases. 

The vertical coordinate for these plots Is normalized In such 
a way that for the typical values of f = 10~‘*sec~^ and Ug = 10 m/sec, 
z may be read directly In meters. For other values of Ug/f , the 
plot must be read as a dimensionless plot and the appropriate 
physical z determined. This shows that for equal Ro the 
neutral boundary layer would become thicker as the equator Is 
approached and/or the geostrophlc wind Is Increased. 

The variances of the three velocity components are shown In 
figures 3 ) and 5 for the same values of Ro . The amplitudes 
of these dimensionless variances decrease as Ro Increases. How- 
ever, if the Ro is Increased by Increasing Ug with Zq and f 
held fixed, then the dimensional value of the variance will 
actually Increase. 

The macroscale A Is shown In figure 6. The solid curves 
give A as predicted by eq. (1^), while the dashed curve gives 
the simplified approximation obtained by holding the slope 9A/3z 
fixed until a maximum value equal to 15 % of the spread of the 
turbulence Is reached. All of the dynamic runs In this report with 
boundary conditions or parameters varying with time have been made 
using eq. (l4). However, the quasi-steady runs shown In figures 1 
though 5 and eqs. (10) and (11) were made with the simpler approxi- 
mation. We do not believe the differences materially affect the 
sensitivity to the parameters shown here. For the specific example 
of the typical diurnal variations, both formulations have been run. 

The results using eq. (l4) are given here and those using the simple 
approximation given in ref. 6. 


Geostrophlc Acceleration 

The next parameter to be Investigated Is to get an idea 

of how important It is to follow the unsteady development of the 
boundary layer rather than confine attention to steady calculations. 
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Figure 2. Influence of Rossby number on the mean wind component 
in the direction normal to the geostrophic wind for the 
conditions of figure 1. 
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Figure 3. Influence of Rossby number on the variance of the 
vertical velocity for the conditions of figure 1. 
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Figure Influence of Rossby number on the variance of the 

velocity in the direction normal to the geostrophic wind for 
the conditions of figure 1. 
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Figure 6. Influence of Rossby number on the turbulent macroscale 
(solid curve as predicted by Eq. (14); dashed curve simplified 
approximation) for the conditions of figure 1. 


To do this, we have made a run where the geostrophlc wind is 
allowed to vary sinusoidally from +10 m/sec to -10 m/sec over a 
period of 24 hours. This yields a maximum value of ~ 6.4 and 
an average value of zero over the 24-hour period. As shown in 
figure 7, there is a significant difference between the steady- 
state values of turbulence intensity and those for the dynamic 
case . 


The characteristic adjustment time in the turbulent correla- 
tion equations, eqs. (4) - (l4), is proportional to A/q . Near 
the surface where A is proportional to z , this characteristic 
time is small — less than a minute for typical conditions at 10 m 
height. Therefore, in the surface laye.r, one can expect the 
turbulent correlations to be in equilibrium with the Instantaneous 
values of the mean flow gradients. However, in eqs. (1) and (2), 
the characteristic time is f which for 45° latitude is of the 
order of 3 hours. Therefore, variations in the geostrophlc wind 
on a time scale of a few hours lead to unsteady effects in the 
velocity gradients through eqs. (1) and (2). This in turn leads 
to variations in the turbulent correlations, as demonstrated in 
figure 7 . 

Near the top of the boundary layer, A/q also approaches f . 
Thus, although a local equilibrium approximation for the turbu- 
lence is quite appropriate near the surface, it is not appropriate 
in the upper levels of the boundary layer when the geostrophlc wind 
is changing on a time scale of a few hours. 

Isolines of constant wind speed in the direction of the geo- 
strophlc wind are shown in figure 8. The Isollnes are nearly 
symmetrical about the 12-hour point where Ug = -10 m/sec. At 
10 m altj^tude, the component of wind shear in this direction is 
0.10 sec \ very close to the steady-state value for the same 
boundary conditions. From figure 1 it may be seen that the wind 
shear may be properly normalized by dividing by f . The dimen- 
sional values quoted in this report are all given for f = 10"‘*sec"^ 

The transverse wind isollnes are shown in figure 9. The 
maximum negative value is shifted to approximately 2 hours after 
the maximum negative value of the geostrophlc wind. In spite of • 
this, the plot tends to confirm the indication of figure 8 that 
the differences in the low-level wind shear between the dynamic 
and steady-state cases are not very significant. The maximum value 
of the vector velocity shear at 10 meters is 0.104 sec~^ for the 
dynamic case in comparison to 0.103 sec"^ for the steady-state case. 
Quite a reduction in the period of the geostrophlc' variation is 
required for this difference to be significant. For the rather 
extreme case of a 12-hour period for a complete cycle, the corres- 
ponding value of the maximum velocity shear remains remarkably 
close, although the height of the boundary layer as marked by the 
height where q^ falls to 10^ of its maximum value is reduced from 
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Figure 9. Isollnes of the component of the mean wind in the direction normal to the 
geostrophlc wind for the conditions of figure 8. 



a value of 2.0 km for the 24-hour period to 1.2 km for the 12-hour 
period case. 


Stability Variations 

Although we expect the dynamics of the changing boundary 
conditions to be Important for most real stability variations In 
the atmosphere, we will first consider the quasi-steady case where 
the surface heat flux Is held constant . Figures 10 and 11 show 
the wind profiles for 4 different stability conditions with Ro 
held constant at 10® . When the bulk Richardson number defined 

by eq. (15) Is positive, the flow Is stable and the turbulence 
damped. As shown In figures 10 and 11, Rl(ig) =0.0135 produces 
a boundary layer which Is significantly smaller than the neutral 
case of RIqqn = 0 , as should be expected for damped turbulence. 
This als'o produces a stronger transverse wind component, as seen 
In figure 11. Although the wind shear at the surface Is reduced 
by the Increase In stability, the wind shear a short distance above 
the surface Is significantly Increased. Thus, from the point of 
view of the wind shear which an aircraft may encounter, an Increase 
In Ri(]_o) quite likely to be associated with Increased wind 
shear. 

As outlined In the previous section, the spec'lflcatlon of 
RIqo) alone Is not sufficient to determine the unstable profile 
(Rx/jn) ^ *^) even In the quasi-steady case. Figures 10 and 11 show 
two^wxnd profiles for nearly the same negative value of Rl(io) 
and two quite different values of the Inversion height . The 
unstable ■ layer , driven by a positive surface heat flux, must be 
capped by a stable, temperature Inversion layer If the boundary 
layer Is to reach a quasi-steady state. The figures show that In 
addition to the shear layer at the surface, a shear layer Is also 
developed In the vicinity of the temperature Inversion. As may be 
noted In the discussion of the boundary conditions In a prior 
section, Zj_ Is not a rigid lid but Is defined as the altitude at 
which the maximum values of 6^ occur as the temperature Inversion 
Is asymptotically approached. In between the surface shear layer 
and the Inversion layer, the layer Is relatively uniformly mixed. 

The vertical velocity variance Is shown In figure 12 for the 
unstable cases but with the normalization which Is appropriate for 
free convection when there Is no mean wind shear. The limiting 
profile for Rl(]_o) ^ , as predicted by our model, was shown 

In reference 8 to agree very well with laboratory simulations of 
this limit. Figure 12 shows that It does not require a very nega- 
tive Ri(io) approach this limiting distribution. The corres- 

ponding values of the Zllltlnkevlch parameter y are also given. 

Next we consider the diurnal variations which may be expected 
to occur on a typical summer day In the mldwestern United States. 
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Figure 11* Influence of Richardson number on the quasi-steady 
distributions of mean wind component In the direction normal 
to the geostrophlc wind for the conditions of figure 10. 
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We have published results for this case (see refs, 5 and 6). In 
figures 13 through 16, the isolines of U , V , , and A are 

shown for a day's cycle for the lowest 1 km altitude. The notation 
for the Isolines is given in Table II. In figures 13b and l4b, we 
also show the vertical profile of the wind components at selected 
times to aid in visualizing the changes. In the plots given pre- 
viously for quasi-steady conditions, it was possible to present re- 
sults in terms of a dimensionless height, zf/Ug . However, in the 
present case, since the dynamics of the dlurnal^varlation is im- 
portant in determining the height of the layer at any given time, 
this nondimenslonallzation is no longer valid. The results of all 
runs involving true time variations in hours will be given in terms 
of altitude in meters. For this run, the geostrophic wind was held 
constant at Ug = 10 m/sec, the upper level temperature gradient 
was held constant at (96/9z)loo = 3°C/km, the surface roughness 
constant at Zq = .01 m , and the surface heat flux was allowed to 
follow an experimentally observed distribution (ref. 19)* Between 
sunset and sunrise, the heat flux is a 'negative constant 
(-0.025°C m/sec) and between approximately 8:00 A.M. and 8:00 P.M., 
the heat flux is positive and nearly follows a sine wave with a 
maximum of 0.25°C m/sec. Over the 24-hour cycle, there is a net 
positive heat flux to the atmosphere which must be removed (in 
simulation of long wave radiation) in order to achieve a periodic 
solution with no long-term heating of the boundary^ layer . In the 
present run, this was accomplished by assuming distributed 

uniformly over the night hours with the following assumed vertical 
distribution : 

• -4 

= -O .703 X 10 cos ('iTz /6000 m)(°C/sec) z <_ 3000 m 

= 0 z > 3000 m 

Figures 13 and l4 show the wind shear to be a maximum in the early 
morning hours when the turbulence, as shown in figure 15 , is a 
minimum. There is a low-level noctural jet which exceeds the 
geostrophic wind by l4^ at a little above 100 m altitude. The 
wind at this time drops from 11.4 m/sec to 3 m/sec at 10 m altitude 
at the same time the principal wind direction is shifting by 
40 degrees. 

A lower scale on the q^ plot (fig. 15) indicates the value 
of the first Zllltlnkevlch parameter as a function of time. 

It is readily observed that the values of q^ occurring when 
Uo = 0 in the morning are quite different from those occurring 
for Uo “ 0 in the evening. 

Figure 16 shows the isolines of the turbulent macroscale A . 
Near the ground, A Increases nearly proportionally to altitude, 
except in the noctural hours when the larger turbulent eddies are 
preferentially damped. This forces the turbulent scale to be less 
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Figure 13a. Isolines of the component of the mean wind in the direction of the 
geostrophic wind for a standard summer day in midwestern United States 
(Ro =10^ , M = Gn = 0). 



Figure 13b. Vertical distributions of u taken at select times 
from figure 13a. 
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Figure l4b. The vertical distributions of v taken at select times from figure l4a. 
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Figure 15. Isollnes of turbulence intensity q‘ 


for the conditions of figure 13 
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Figure 16. Isolines of turbulence macroscale A for the conditions of figure 13. 


at the same altitude than It Is during the daytime hours and Is 
reflected In figure l6 as a rise In the Isollnes during the noctur- 
nal hours. At moderate altitudes, the scale remains relatively 
large even during the early morning hours. However, as noted from 
figure 15, this Is associated with very low Intensity turbulence. 

So the flow Is very nearly laminar even though the scale of what 
turbulence remains Is large. 

Very little difference Is observed If only the upper-level 
temperature gradient Is changed. A run was made with conditions 
Identical to the preceding run with the exception of (90/3z)|^ 
changed from 3°C/km to the Isothermal value of 9.87°C/km. The 
resulting wind and turbulence distributions showed no significant 
differences from those given In figures 13 through l6. Even the 
value of the height of the upper-level Inversion which was achieved 
when the flow became periodic was surprisingly close to that of the 
previous run. 

Figures 17 through 20 show the variation caused In the diurnal 
fluctuations If the afternoon positive heat flux Is reduced In 
addition to Increasing the upper-level lapse rate. For this run, 
the surface heat flux was taken as exactly symmetrical In time so 
there Is no need for a radiation flux divergence term to make the 
boundary layer periodic. We specified w0|o = .025 sln(Tr/12 ) (t-8 
hrs) and Zq = 0.1 m. Also the upper-level potential temperature 
gradient Is Increased from 3°C/km to the Isothermal value of 
9.87°C/km. Other conditions are the same as for figures 13 through 
16. In both cases, the flow Is permitted to run through three 
complete cycles to eliminate any Influence of the Initial condi- 
tions. The third cycle appears quite close to periodic and Is the 
result shown. 

There Is considerably less time variation In the turbulence 
level shown In figure 19. There Is neither sufficient heating In 
the afternoon to force the large ballooning of the turbulence 
Intensity to higher levels seen In figure 15 nor sufficient cooling 
at night to strongly damp the turbulence. However, as seen In 
figures 17 and l8, the wind distribution Is markedly different from 
either the standard day variation of figures 13 and l4 or the 
neutral, steady-state profiles of figures 1 and 2. The low-level 
nocturnal jet Is actually Increased to 13 -1 m/sec and occurs at a 
somewhat higher altitude, as seen by comparing figures 13 and 17. 
This Is an example of how difficult It Is to track through all the 
nonllnearltles In the equations and anticipate the results of some 
particular change. 

The macroscale In figure 20 Is reduced throughout the day to 
a value which Is even less than that for a steady-state, neutral 
boundary layer. This may be a result of the Imposed temperature 
Inversion which was Initially started at 100 m altitude. In the 
standard day run (figs. 13-16), the temperature Inversion was 
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Figure l8b. 


Vertical distributions of v 


taken at select times from figure l8a. 





raised by the turbulence to approximately 3 km. The reduced 
surface heat flux of the present case only pushes the Inversion 
to approximately 500 m altitude. This is much less than the height 
of the neutral, steady-state boundary layers shown in figures 11 
through l6. It appears to permit the natural inertial oscillations 
to play a more Important role, particularly in the transverse wind 
component V . 


Influence of Barocllniclty 

The approximation of horizontal, spatial Inhomogeneity is 
seldom completely realized In the atmosphere. One effect which 
perhaps Is most often encountered Is that of horizontal temperature 
gradients. Fortunately, this effect can be Incorporated within the 
one-dimensional model by permitting the pressure gradient, repre- 
sented In eqs. (1) and (2) in terms of geostrophlc velocities, to vary 
with altitude. This "thermal wind" variation In the geostrophlc 
wind can be obtained directly from the geostrophlc balance by 
differentiating and using the equation of state 
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with X in the direction of Ug at the top of the boundary layer. 
If temperature advection Is to Be Included In the energy equation, 
then the source term in eq. (3) should Include a term of the form 
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( 22 ) 


We will consider the direct effect of the geostrophlc wind 
shear In the momentum equation separately from the effect of temp- 
erature advection. Figures 21 through 2^ show the results for 
M = 1.25 . To Include both the dynamics and the Influence of 
different values of Bq , we have allowed 3q to vary from 0 to 
2 tt and back to 0 In a period of 12 hours. The magnitude corres- 
ponds to a horizontal temperature change of 3°C In 100 km which 
appears as a not uncommon value in nature. The time variation 
corresponds to a relatively fast moving temperature anomaly, such 
as a thunderstorm at a distance. Again, the flow was permitted to 
become relatively periodic by allowing 3 full periods of variation 
before the results shown In the figures were recorded. For this 
first set of figures, the temperature Is uniform between the surface 
and the inversion layer at approximately 1250 m altitude. 
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Figure 21. Isollnes of u for a thermal wind varying periodically In direction 
(Ro = 10® , 96/3zL = 3°C/km, M = 1.25 , = 0 , Rinn^ = 0)- 
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The maximum U velocity, 13*5 m/sec, occurs approximately 
3 hours after the thermal wind Is directed in the direction of the 
geostrophlc wind (3 q = 0)* But 3 hours after the thermal wind is 
next directed in the same direction (Bo = Sir), U is near a 
minimum. Actually, from eqs . (1) and (2) we see that the U 
velocity responds to variations in Vg and, since 3Vg/9z follows 
sin Bt , this contributes a positive pressure gradient in eq. (1) 
between 0 < Bq ^ ^ negative pressure gradient between 

7T < Bq < 2it . Thus it is consistent for a maximum to occur in U 
when Bq passes through tt as it is Increased and for a minimum 
to occur in U when Bq passes through tt as it is decreased. 
This particular time variation of driving conditions is not 
expected to occur long enough in nature to make the real flow 
periodic, but it does demonstrate how Important the time history 
of the flow is in determining the response of the flow to instan- 
taneous values of the driving conditions. 

The maximum value of V , 7 m/sec, occurs near the 10-hour 
mark after the combination of a relative minimum in U and a 
maximum in Ug has produced a favorable pressure gradient (see 
eq. (2)) for a period of time. There is a region of slightly 
negative V (approx. -1 m/sec) between 1000 m and 1250 m altitude 
which does not show in figure 22. The turbulence Isollnes in 
figure 23 show no strong variation over the period. Neither do 
the Isolines for A in figure 24, although there is a slight maxi- 
mum occurring in A near the time that V has its local maximum. 

The maximum wind shear at or above 10 m altitude computed by 
taking the square root of the sum of the squares of 3U/9z and 
3V/3z gives a value of 0.l4 sec~^ , somewhat less than the maximum 
value (0.l8 sec'M that occurs in the typical day calculations of 
figures 13 and 14. However, in this case the shear occurs over a 
somewhat wider altitude range. Perhaps the more Important differ- 
ence is that in the earlier figures the strong local stability 
caused by surface cooling was the dominant mechanism while, in the 
present case, there is a neutral, vertical temperature variation 
in the region of maximum shear. This suggests two different ways 
in which the wind shear may be Increased: either by increasing 

the horizontal wind gradients to increase the driving thermal winds 
or by applying a vertical temperature gradient which Increases the 
stability of the flow under the same driving pressure gradients. 

Figures 25 through 28 show the results of taking the same 
conditions as those for figures 21 through 24 and changing only the 
upper-level lapse rate from 3°C/km to 6°C/km. As seen by comparing 
figures 21 and 25, the maximum U velocity is Increased to l6.6 
m/sec from 13.6 m/sec. Also from figures 22 and 26, the maximum V 
is increased from 7 to 10 m/sec. However, the resulting maximum 
shear still occurs near the surface with the value at 10 m altitude 
unchanged at 0.l4 sec"^. 
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Figure 28. Isollnes of turbulence macroscale A for the conditions of figure 25 




If, Instead of increasing the stability, the magnitude of the 
driving thermal wind is increased an order of magnitude corres- 
ponding to 3°C in 10 km (M = 10) but other conditions remain the 
same as in figures 21 through 24, the maximum velocity in the U 
direction is increased to 25.4 m/sec and that in the transverse 
direction to 24.6 m/sec while the turbulence increases to 

12.5 (m/sec) The maximum wind shear for this case in increased 
to 0.33 sec~^. The maximum velocities associated with this case 
are close to the maximum gust velocities typically associated with 
the passage of thunderstorms. It might be anticipated that the 
maximum velocities and also the wind shear could be further 
Increased by applying stable, vertical temperature gradients to 
this case. However, the large levels of turbulence associated 
with this case tend to rapidly mix any vertical temperature grad- 
ients before they significantly damp the turbulence. Thus, only 
quite modest further Increases in maximum velocities and maximum 
wind shear are achieved by imposing any reasonable vertical 
thermal boundary conditions on this problem of very large M . 

The Influence of temperature advectlon is closely associated 
with the stabilizing effect of the vertical thermal boundary layer 
conditions. The computer run associated with figures 25 through 
28 was rerun with temperature advectlon added. No other changes 
were made in the input conditions. The major difference in the 
computed results is that the advectlon destroys the periodicity of 
the output. Due to the nonlinearity in the source term, eq. (22), 
there is a net energy sink. So there is continued cooling in the 
mixed layer between the surface and the inversion layer. Figures 
29a and 29b show the difference in the wind and temperature mean 
profiles. two complete periods after the run is started with the 
same initial conditions. 12 hours later, run 16 would again 
predict the same profiles, but run 17 would have a further 
decreased mixed layer temperature and an even stronger shear layer 
in the neighborhood of 300 m altitude. 

The results of our runs show that the Influence of barocllnl- 
clty is strongly affected by the time history of its development 
on the time scale over which it is likely to occur in nature. 

Thus, in any attempt to simulate conditions in the barocllnlc 
atmospheric boundary layer at a given time and place, it is neces- 
sary to trace the development of the boundary layer for some time 
preceding the desired time. In the case of our attempt (ref. 9) 
to simulate the conditions at JFK airport at the time of the 
Eastern Airlines crash in June 1975, we started the simulation 
l6 hours in advance of the desired time. The detailed distribution 
predicted and presented in reference 9 will not be repeated here. 
The general conclusion was that a wind shift of the order of 
15 m/sec within a 100 m altitude change was consistent with the 
meteorological conditions leading up to and at the time of the 
crash. Shear of this same order of magnitude is given in several 
figures of this report, but the wind shift over any 100 m altitude 
change is less. 
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Figure 29a. Influence of temperature advectlon on the profiles of 
figures 25 and 26 at a set time, 2H hours, after initialization 
of the run (Run 16 without advection and Run 17 with advection) 
(a) mean wind components. 
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Figure 29 b. Influence of temperature advectlon on the profiles of 
figures 25 and 26 at a set time, 2k hours, after initialization 
of the run (Run I 6 without advection and Run 17 with advectlon); 
(b) potential temperature distribution. 
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Radiation Flux Divergence 

The distributions of wind and turbulence may be strongly 
affected by the direct energy source or sink within the boundary 
layer by thermal radiation flux divergence. This has already been 
Indirectly demonstrated by the influence of the upper-level temper- 
ature gradient on some of the previous solutions, since radiation 
is important in maintaining this stable upper-level temperature 
gradient. Here we would like to carry this a step further In 
demonstrating that under some conditions radiation flux divergence 
can govern the diurnal cycle. 

In places where the earth’s surface has a very high heat 
capacity, the surface temperature may be taken as essentially 
constant over time periods of the order of a day. This Is the case 
for the ocean. Under the conditions of a constant surface tempera- 
ture, any thermal diurnal variations will be driven by the direct 
energy absorption or emission from the gas In the boundary layer. 
Figures 30 through 33 simulate such a possible diurnal variation. 
For this case, the geostrophlc wind was maintained constant at a 
typical value of 10 m/sec. The upper-level potential temperature 
gradient was maintained at 3°C/km, the aerodynamic surface rough- 
ness was set at 0.10 m, and the surface temperature held fixed. 

A radiation source term of the form 


= .000556 sin (^) cos ( ( e"^/500 ^2^c>c/sec ) (23) 

was assumed. This is the type of radiation flux divergence which 
might be expected In a boundary layer with a relatively strong 
decrease In water vapor with altitude. The periodic variation with 
time was again chosen to eliminate initial condition dependency. 

The resulting gravitational instability in the vicinity of the 
top of the heated layer forces this boundary layer to be relatively 
thick (approx. 1.8 km). We have only shown the lowest 1 km In 
figures 30 through 33. The turbulence Intensity plot (fig. 32) is 
quite similar to that for the typical mldwestern summer day (fig. 
15). Although the stability in the surface layer Is 12 hours out 
of phase with respect to the earlier run (l.e., the night time 
surface layer is unstable while the day time surface layer is 
stable), the maximum turbulence Intensity Is dominated more by the 
stability In the upper regions of the boundary layer which follows 
the same time dependence as that for figure 15. This upper region 
stability also dominates the turbulent scale variation In figure 33. 

The maximum U wind component occurs In the afternoon, while 
the maximum V component occurs In the late morning. The maximum 
wind shear above 10 m altitude Is 0.12 sec“^, occurring at approxi- 
mately 15 hours. 
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CONCLUDING REMARKS 


We have not attempted to predict results for all possible 
combinations of the governing parameters. In fact, we have shown 
that even if this were done it would not exhaust the possibilities 
for wind and turbulence profiles which might exist, due to the 
importance of the time history of the governing parameters. Even 
within the limitations of our one-dimensional approximation, there 
is Indeed a rich assortment of possible wind and turbulence pro- 
files. What we have shown here is a sampling of the possibilities 
with some Indication of the Influence of the governing parameters. 

A partial summary of the turbulence intensity for different 
conditions is shown in figure 3^. Refer to Table III for the key 
to the different conditions. The turbulence is normalized by the 
maximum wind velocity between the ground and 1 km and plotted as 
a function of stability. As indicated in the previous section, it 
is not possible to reduce the detailed Influences of stability to 
a single parameter. But for purposes of indicating only the rough 
dependence of stability, we have chosen a single parameter which 
does partially combine the influence of the surface heat flux 
which controls stability in the surface layer and the Influence of 
the upper-level potential temperature gradient which Influences 
the elevation of the inversion layer which caps the boundary layer. 
This parameter is a bulk Richardson number based on the velocity 
and potential temperature differences between the surface and 1 km. 
The hysteresis loops presented for the different periodic cycles 
considered indicate the type of scatter to be expected in at- 
tempting to predict the fluctuations in this manner. As the bulk 
Richardson number approaches -» , the conditions for free con- 
vection, the normalized q will approach “ . On the stable side. 
Increasing f^l(iooo) tends to decrease the normalized q whether 
the increase in stability is caused by Increased surface codling or 
a decrease in the upper-level inversion height. 

A similar rough summary of the maximum wind shear occurring at 
or above 10 m altitude in the boundary layer is given in figure 35- 
There is a fairly consistent trend of increasing shear with in- 
creasing stability, although there is a relatively strong Influence 
of thermal wind which tends to mask this trend to some extent. 

The normalization of q by the maximum velocity within the layer 
tends to roughly compensate for the thermal wind influence in 
figure 3^ . Of course, the wind shear in a profile at a given time 
cannot be characterized by a single number. The actual changes in 
airspeed which an aircraft will encounter will also depend on the 
vertical depth over which this large shear exists. This must be 
determined by looking at the detailed profiles. Figure 35 does 
indicate the rough rule of thumb that as either horizontal temper- 
ature gradients Increase or vertical stable temperature gradients 
increase, problems with wind shear may be expected to increase. 
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Figure 3^ • Summary plot of normalized total velocity fluctu 
ation as a function of the bulk Richardson number between 
the surface and 1 km altitude. 
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Figure 35* Summary plot of maximum wind shear at or above 10 m 
altitude as a function of the bulk Richardson number between 
the surface and 1 km altitude. 
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In general, the wind shear decreases with altitude in the 
lower portion of the boundary layer. This Is shown in figure 36 
to supplement the summary Information In figure 35. The solid 
curve shows the distribution of wind shear with altitude for the 
neutral, barotroplc distribution of figure 1. The bars represent 
the change in wind shear caused by two orders of magnitude change 
in Ro. The dashed curves represent the bounds in the variation of 
wind shear caused by expected values of Ri. In general, increas- 
ing R1 will tend to increase wind shear, as seen in the previous 
figure, but the minimum wind shear need not coincide with the 
minimum Ri due to hysteresis effects. Barocllnicity also will, in 
general, increase wind shear. It is not possible to correlate the 
variations indicated in figure 36 with time of day. Nocturnal 
hours do tend to be more stable (Ri > 0) and afternoon hours more 
unstable (Ri < 0), but the barocllnicity is associated with 
changing weather and does not correlate with time of day. A 
companion curve like figure 36 for the turbulence Intensity as a 
function of altitude would be too misleading, since there is no 
consistent trend of turbulence Intensity with altitude over para- 
meter space. Under stable conditions, the turbulence will 
decrease, as seen in the early morning hours of figure 15, with 
altitude, but for unstable conditions at low levels, it will 
Increase strongly with altitude, as seen in the afternoon hours of 
the same figure. 

Many cases in nature where wind shear and turbulence are a 
problem to aviation are further complicated by the added dimension- 
ality of the phenomena. For example, a strong downdraft may be 
associated with flow from a thunderstorm. To Include the effect of 
such vertical velocities, it is necessary to at least extend our 
problem to a two-dimensional, unsteady program. This has been done 
during the time period of this contract under separate ERA and Navy 
contracts. This opens up the possibility of predicting the wind 
shear associated with strong local cold fronts, but this is beyond 
the scope of the present Investigation. 

We believe there are two tasks remaining to be accomplished 
in our effort to determine the feasibility of predicting the general 
meteorological conditions which lead to wind shear and turbulence 
conditions most hazardous for aircraft takeoffs and landings. One 
is a detailed investigation of the ensemble of aircraft accidents 
identified in Part 2 of this report and, second, a detailed program 
of comparison with turbulence and wind profiles recorded by airport 
instrumentation for the purpose of determining the accuracy of our 
model predictions. This combination of looking at cases which have 
proven to be a problem in the past, together with a careful deline- 
ation of the error Involved in our predictions by additional 
comparisons with field observations, will allow us to determine the 
feasibility of predicting dangerous situations in advance. We 
believe that an interaction between our modeling efforts and the 
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meters 



Figure 36. Summary plot of wind shear as a function of altitude showing the influence 
of the expected range of variations in Ro, Ri and M. 


measurements program also would be of benefit to the measurements 
program by Indicating possible modes for organizing the data. 
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TABLE II 


NOTATION FOR ISOLINE PLOTS 
(For Figs. 8-9 > 13-28, and 30-33) 


Contour § 

represents 

v/v 

g 

q 

^(m^/sec^) 

A(m) 

0 

-0.9 

-0.2 


0 

1.0 

1 

-0.6 

-0.1 


0.1 

10.0 

2 

-0.3 

0.05 


1.0 

20.0 

3 

+0.1 

0.1 


1.5 

50.0 

4 

0.3 

0.2 


2.0 

75.0 

5 

0.6 

0.3 


2.5 

100. 

6 

0.9 

0.4 


3.0 

200. 

7 

1.2 

0.6 


4.0 

300. 

8 

1.5 

0.9 



400. 

9 

■< 95 % 

of maximum 

value 
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TABLE III 

KEY FOR SUMMARY PLOTS (FIGS. 3^ AND 35) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


For the conditions corresponding to figures 13 - l6 
For the conditions corresponding to figures 17 - 20 


For the conditions corresponding to figures 30 - 33 
For the conditions corresponding to figures 21 - 24 


For the conditions corresponding to figures 21 - 24 but 
Mq increased from 1.25 to 10 

For the conditions corresponding to figures 25 - 28 


Neutral steady state 
Neutral steady state 
Neutral steady state 
Stable atmosphere 
Stable atmosphere 
Stable atmosphere 
Unstable atmosphere 
Unstable atmosphere 


(Ro 


10 

5 

R1 

= 

0) 

(Ro 

= 

10 

6 

3 

R1 

= 

0) 

(Ro 

= 

10 

7 

3 

R1 

= 

0) 

(Ro 

= 

10 

6 

3 

R1 

(10 

) 

(Ro 

— 

10 

vJ 

3 

7 

R1 

(10 

) 

(Ro 

~ 

10 

l 

3 

R1 

(10 

) 

(Ro 

= 

10 

U 

3 

7 

R1 

(10 

) 

(Ro 

= 

10 

( 

3 

R1 

(10 

) 


0.014) 

0 . 0036 ) 

0.047) 

0.0128, z^f/u = 0.21 
0.035, z^f/u = 0.16 


For conditions simulating JFK weather at the time of EA Flight 
66 crash (fig. 15, ref. 9) 
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WIND SHEAR AND TURBULENCE AROUND AIRPORTS 


PART 2. SELECTION OF AN ENSEMBLE OF AIRCRAFT ACCIDENTS 
FOR WHICH WIND SHEAR OR TURBULENCE APPEARED TO 
BE A CONTRIBUTING FACTOR 


W. S. Lewellen and Guy G. Williamson 


INTRODUCTION 


This part of the study is an investigation of a set of 
accidents where wind shear and/or turbulence were significant 
factors. The study was conducted in order to select cases which 
can be used to evaluate the feasibility of predicting when these 
weather conditions will be hazardous to aircraft operation. 

It is well known that wind shear and gusts Increase the 
difficulty of aircraft operation, particularly during landing and 
takeoff (ref. 20). Gusts or rapid changes in the wind vector as 
a function of time produce undesirable aircraft accelerations 
which must be compensated for by control Inputs, thus increasing 
the pilot's workload. Wind shear or changes in the wind vector 
as a function of spatial location produce similar results. If the 
change . in wind vector is small or occurs slowly, it does not 
significantly affect aircraft operation. However, rapid large 
changes can cause crashes, regardless of what the pilot and/or 
automatic controls do. It is therefore desirable to know 

a) what magnitude of wind shear and turbulence pose threats 
to aircraft safety, and 

b) when and where such threats exist. 

Answers to the first part of the question have been based on past 
operating records and simulator studies. The solution to the 
second part requires direct measurement of these variables or 
reliable predictions of their likely existence. 

Measurement of these variables throughout the landing and 
takeoff space requires elaborate, expensive equipment and it is not 
considered feasible for many airports at the present time. On the 
other hand, the reliability of predicting these conditions must 
still be shown. 

One step in evaluating a predictive technique is to compare 
its results with known weather conditions which probably caused a 
real aircraft accident. If the technique can predict weather 
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conditions similar to those at the time of the accident, then it 
can conceivably be used to warn aircraft operations personnel of 
impending unsafe flight conditions. In order to perform this 
evaluation, it is necessary to select a number of real accident 
cases for use in the comparison. 

Turbulence and/or wind shear are associated with several 
different weather conditions. An obvious hazardous weather condi- 
tion is in the vicinity of thunderstorms. Here one expects strong 
turbulence, horizontal wind shear, and downdrafts. A not so 
obvious condition is when very stable temperature gradients prevail 
and high shears and low turbulence are possible. This type of 
condition can also have fog associated with it. The combination of 
low visibility due to fog and unexpected wind shear can give rise 
to aircraft accidents in which the wind shear factor goes unnoticed. 
These stable conditions can also occur during the passage of a cold 
front with the cold air slipping under the existing warm air mass. 

A third condition of Interest occurs when a strong wind condition is 
combined with a sizable change in surface roughness. This combina- 
tion, as exemplified by a sea breeze passing over a hilly Island, 
can'produce hazardous gusty conditions. 

The selection of accidents for future study has been conducted 
with an eye towards obtaining cases representative of a number of 
different conditions. Such a study would permit a catalogue of 
hazardous weather conditions to be constructed with the complete 
model description of the detailed profiles of wind magnitude and 
direction and ensemble-averaged, turbulent fluctuations for each 
case. The selection process is further described and the resulting 
cases listed in the following section. 
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SELECTING ACCIDENT CASES 


The starting point for the selection of accident cases was 
the National Transportation Safety Board (NTSB) aircraft accident 
data stored on magnetic tape. A description of each accident Is 
generated and stored In coded form as specified in reference 21. 

The important characteristics of the accidents for this study are 

1. Operational phase of flight 
(landing, takeoff, cllmbout, etc.) 

2. Type of accident 

(stall, spin, hard landing, collision, etc.) 

3. Probable causes of accident 
(weather, pilot, terrain, etc.) 

Since the study Is concerned with low altitude wind shear and turb- 
ulence, the operational phases were limited to takeoff, landing, 
and missed approach. Similarly, the types of accidents which are 
likely to be caused by these weather conditions are hard landing, 
collision with ground obstacles, overshoot, undershoot, and stall. 
Finally, the probable causes or factors Involved In the accidents 
are particular weather conditions and pilot error. Accidents due 
to certain pilot errors have been Included because (a) pilot errors 
themselves could be Induced by wind shear and turbulence, or (b) 
the effects of wind shear could be mistaken for pilot error. The 
above-mentioned accident characteristics are summarized In Table 
IV. Note that "wind shear" Is not a codeable cause/factor, so 
other possibly related cause factors must be used. 

The NTSB aircraft accident data bank was queried for the time 
period 1970-1974. The search was done in several parts with some 
accidents appearing In more than one part. Adjusting to approxi- 
mately remove the redundancy from the accident count results In 
1160 accidents In which wind shear or turbulence may be Involved. 
Obviously, this large number of accidents should be further 
screened In order to obtain a more select set for further study. 

The next cut was based on the following observations: 

1. The stronger the wind shear or turbulence, the greater 
the aircraft upset and, therefore, seriousness of the 
accident . 

2. The accident data files for commercial accidents are 
more comprehensive than those for noncommercial ones. 

Therefore, noncommercial accidents listed as less than serious 
were, with a few exceptions, deleted from the candidate list. The 
exceptions are cases where either the terrain or weather condi- 
tions made them particularly interesting. This left 47 takeoffs 
and 85 landings, a total of 132 candidate accidents. 
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TABLE IV. 


Cause/Factors , Type of Accident, and Operational Phases 
where Wind Shear and/or Turbulence Could Have Contributed 

to the Aircraft Accident 


Cause/Factors 


Sudden wind shifts 
Unfavorable wind conditions 
Downdrafts, updrafts 

Improper compensation for wind conditions 
Fog 

Thunderstorm activity 

Clear air turbulence 

Palled to maintain flying speed 

Misjudged distance, speed, altitude 

Delayed initiating go-around 

Palled to Initiate go-around 

Continued flight into known area of severe turbulence 
Lost /disoriented 
Misused flaps 

Spontaneous Improper action 

Palled to maintain directional control 

Vortex turbulence 


Type of Accident 
Hard landing 

Collision with ground - controlled 
Collision with ground - 
uncontrolled 

Collision with obstructions 

Overshoot 

Undershoot 

Stall 

Stall - spin 
Stall - spiral 
Stall - mush 
Turbulence 

Uncontrolled altitude deviation 


Operational Phase 

Landing - final VPR 

Landing - final IFR 

Landing - initial 

Landing - level off/touchdown 

Landing - go-around 

Landing - missed approach 

Takeoff - initial climb 
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The next cut from this list was made on the basis of the 
complexity of modeling the local weather conditions. In particular, 
noncommercial accidents occurring In areas where there are gross 
changes In the features of the terrain or occurring within thunder- 
storms were excluded due to modeling difficulties. This left 53 
accidents to be studied In more detail. 

The files for 40 of these accidents were reviewed at the 
Office of Public Inquiries In Washington, and l8 cases were selected 
for further study. The others were excluded primarily due to 
complex terrain and crosswind effects as opposed to wind shear and 
turbulence. In addition, accidents Involving pilots with less than 
1000 total hours and less than 100 hours In that type of aircraft 
were excluded In order to reduce the probability that pilot Inex- 
perience was a cause. It Is anticipated that a few of the remain- 
ing accidents will be eliminated on the basis of further review. 

For accidents occurring to scheduled flights, the file can run 
several hundred pages and usually contains all the data necessary 
for local weather analysis. However, the vast majority of the 
accidents occur on nonscheduled flights and the files associated 
with these are usually Just a few pages. For further study of 
these cases. It will be necessary to obtain the appropriate 
weather data from the National Weather Service. Terrain data In 
the vicinity of some of the airports will also be required. The 
particulars for the final set of accidents selected are given In 
Table V. 
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TABLE V. ACCIDENTS FOR FURTHER INVESTIGATION 


File No. 

Location 

Date 

A/C No. 

Commercial 
or Noncomm. 
Flight 

Accident 

Type 

Operational 

Phase 

Cause/ 

Factor' 

A-0004 

MA . , Boston 

12 / 17/73 

EC-CBN 

C 

Hit obstruction 

Landing 

S.W. 


NY-JFK 

6/24/75 

N8845E 

C 

Hard landing 

Landing 

D-U 

1-0041 

Missouri , 
St .Louis 

7 / 23/73 

N4215 

C 

Hard landing 

Landing 

TSTM, 

U.W. 

1-0028 

Tennessee , 
Chatanooga 

11 / 27/63 

N3323L 

C 

Undershoot 

Landing 

U.W. 

1-0047 

NY-JFK 

12 / 12/72 

N788TW 

C 

Hit obstruction 

Landing 

S.W. 

3-0001 

NY-LaG. 

1/4/71 

N7 

C 

Hit obstruction 

Landing 

U.W. 

3-0001 

Penna . , 
Johnstown 

1 / 6/74 

N125AE 

C 

Stall 

Landing 

D-U 

3-4528 

Puerto Rico, 
Culebra 

12 / 21/71 

N589JA 

c 

Hit ground 

Landing 

U-D 

3-0617 

N. Carolina, 
Morrlsvllle 

4 / 2/70 

N401RA 

NC 

Hit trees 

Landing 

S.W. 

3-3086 

Conn . , 
Madison 

8 / 26/74 

N5065B 

C 

Undershoot 

Landing 

D-U 

3-1898 

Washington, 

Mead 

5 / 19/74 

N55 

NC 

Undershoot 

Landing 

D-U 


Standard cause/factor terms: S.W. = sudden wind shift; D-U = downdraft-updraft; 

U.W. * unfavorable wind; TSTM = thunderstorm 



Pile No. 

Location 

Date 

A/C No. 

Commercial 
or Noncomm. 
Flight 

Accident 

Type 

Operational 

Phase 

Cause/ 

Factor 

3-0780 

Georgia, 
St .Marys 

3/26/70 

N7332T 

C 

Stall 

Takeoff 

U.W. 

1-0050 

Virgin Is., 
St . Thomas 

12/10/70 

N3417 

' C 

Hard landing 

Landing 

U.W. 

1-0012 

Calif. , 

Los Angeles 

1/16/74 

N757TW 

C 

Hard landing 

Landing 

Fog 

1-0005 

Kansas, 

Wichita 

3/3/74 

N2307 

C 

Overshoot 

Landing 

Rain 

1-0041 

Alaska, 

Sltkinak 

10/2/70 

N9117 

C 

Stall 

Takeoff 

U.W. 

1-0002 

Calif. , 
Ontario 

3/31/71 

N3166 

C 

Hit ground 

Landing 

Fog 

1-0008 

Alaska, 
N. Juneau 

9/4/71 

N2969G 

C 

Hit ground 

Landing 

Fog 


Note: The files of those accidents listed above the dashed line were reviewed in some 
detail while those below the dashed line were selected on the basis of the 
accident brief only. 
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